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A new calculation of the p/p ratio in cosmic rays is compared to the recent PAMELA data. 
The good match up to 100 GeV allows to set constraints on exotic contributions from thermal 
WIMP dark matter candidates. We derive stringent limits on possible enhancements of the WIMP 
p flux: a mwiMP = 100 GeV (1 TeV) signal cannot be increased by more than a factor 6 (40) without 
overrunning PAMELA data. Annihilation through the W + W~ channel is also inspected and cross- 
checked with e + /(e _ + e + ) data. This scenario is strongly disfavored as it fails to simultaneously 
reproduce positron and antiproton measurements. 



PACS numbers: 95.35.+d,98.35.Gi,11.30.Pb,95.30.Cq 

The cosmic ray (CR) antiproton and positron fluxes 
are considered as prime targets for indirect detection 
of galactic dark matter (DM). A deviation from the 
predicted astrophysical background has been searched 
for mostly at low energy (e.g., [l|). However, in 
some scenario, heavy Weakly Interacting Massive Par- 
ticle (WIMP) candidates — either from annihilation of 
the lightest supcrsymmetric species, or from the light- 
est Kaluza-Klein particles in universal extra dimen- 
sions — should be able to provide sizeable fluxes beyond 
a few GeV Q . The PAMELA collaboration has recently 
published the cosmic ray antiproton to proton ratio in 
the hitherto unexplored ~ 1 — 100 GeV energy range 
We find that the background flux, yielded by standard 
astrophysical processes, can explain the data up to high 
energies. Adding a contribution from the annihilation 
of a generic WIMP dark matter halo, we derive strin- 
gent upper limits on possible boost factors of the exotic 
p flux. In addition, in some scenarios, the WIMPs mostly 
annihilate into W + W~ pairs, hence giving rise to a copi- 
ous amount of hard positrons. Recent calculations, in the 
light of the new PAMELA measurement of the positronic 
fraction in CRs, showed that the secondary production 
alone could marginally explain the data Should a 
heavy WIMP be required to better match the positron 
flux, we also check the viability of such models through 



the combined constraints from the e + and p fluxes as also 
done in Q . The secondary p flux provided in this paper 
is an improved calculation of that presented in [19| , to 
which we refer for a thorough discussion of the ingredi- 
ents and the technical details. Antiprotons arc yielded 
by the spallation of cosmic ray proton and helium nuclei 
over the interstellar medium, the contribution of heav- 
ier nuclei being negligible [2fj|. Even if only p and He 
with kinetic energy larger than 6 m p can produce p, a 
good description of the p and He interstellar (IS) fluxes 
is mandatory to correctly provide the p/p ratio in the 
0.1 GcV-100 GeV range. Following 0|, we model the 
proton and helium IS fluxes as 



$ = A/3 n R 



Pi R-P2 



- 1 (GcV/n)- 



(1) 



where R is the rigidity of the particle. The param- 
eterization for the fluxes below 20 GeV/n are taken 
from the reanalysis of the 1997 to 2002 BESS data 
ill: {A,P U P 2 } = {19400,0.7,2.76} for H and 
{7100, 0.5, 2.78} for He. For the high energy range, 
the combined fit of AMS-01 0, El Q, BESS98 
and BESS-TeV Q demodulated data respec- 
tively give {A, Pi, P 2 } = {24132, 0., 2.84} for H and 



0., 2.85} for He. The two fits connect smoothly 
at 20 GeV/n. Compared to 19(, we also improve the 
calculation of the tertiary mechanism 27] . The Ander- 
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FIG. 1: Left panel: IS antiproton flux for the B/C best fit model and two parameterizations of the production cross section. 
Right panel: propagation uncertainty envelopes of the IS p/p ratio for the same production cross sections as in the left panel. 
All data are demodulated using the force-field approximation: AMS 98 jfj, IMAX 92 @], CAPRICE 94 || .WIZARD-MASS 91 
[3, CAPRICE 98 BESS93 ELbESS 95+97 [H J3ESS 98 0, BESS 99 and 2000 Q, BESS 2002 BESS Polar 0, 
WIZARD-MASS 1[13], HEAT-JHa], and PAMELA 0]. 



son prescription [28| is used, as described in As 
a result, the low energy tail is more replenished, lead- 
ing to a larger flux. The framework used to calculate 
cosmic ray fluxes is the diffusion model with convec- 
tion and reacceleration. The transport parameters are 
fixed from the boron-to-carbon (B/C) analysis (HB] and 
correspond to i) the diffusion halo of the Galaxy L; ii) 
the normalization of the diffusion coefficient Kq and its 



slope 5 (K(E) = K f3R s 



the velocity of the con- 



stant wind directed perpendicular to the galactic disk 
V c = ±V c e z ] and iv) the reacceleration strength medi- 
ated via the Alfvenic speed V a . Strong degeneracies are 
observed among the allowed parameter sets |30|. but it 
has a limited impact on the secondary p flux [l9(. At 
variance, the corresponding DM-induced p flux suffers 
large propagation uncertainties [3l|; this is also the case 
for the secondary and primary positron fluxes 0, [j^ ]. 
Throughout the paper, the fluxes will be shown for the 
B/C best fit propagation parameters, i.e. L = 4. kpc, 
A'o = 0.0112 kp^Myr" 1 , S = 0.7, V c = 12. km s" 1 and 
V a = 52.9 km s" 1 0. 

The secondary IS p flux is displayed in the left panel 
of Fig. [T] along with the data demodulated according to 
the force-field prescription. We either use the DTUNUC 
[T9I ] p production cross sections (solid line) or those dis- 
cussed in 0, 20 1 (dashed line). The differences between 
the two curves illustrate the uncertainty related to the 
production cross sections, as emphasized in [l9j |. where 
a careful and conservative analysis within the DTUNUC 
simulation settled a nuclear uncertainty of ~ 25% over 
the energy range 0.1 — 100 GeV. The conclusion is sim- 



ilar here, although the two sets of cross sections differ 
mostly at low energy. In the right panel, along with the 
demodulated p/p data, we show the curves bounding the 
propagation uncertainty on the p calculation based cither 
on the DTUNUC [19( p production cross sections (solid 
lines) or those borrowed from 0] (dashed lines) . The un- 
certainty arising from propagation is comparable to the 
nuclear one [l9| . From a bare eye inspection, it is evident 
that the secondary contribution alone explains PAMELA 
data on the whole energetic range. It is not necessary to 
invoke an additional component to the standard astro- 
physical one. 

Motivated by the accuracy of our predictions and their 
well understood theoretical uncertainties, as well as by 
the good statistical significance of PAMELA data, we 
derive limits on a possible exotic component. We fo- 
cus on the high energy part of the p/p ratio, where so- 
lar modulation does not play any role We assume 
an additional component of antiprotons produced by an- 
nihilation of WIMPs filling the dark halo of the Milky 
Way. Their distribution is taken as a cored-isothermal 
sphere with local density pq = 0.3 GeV cm~ 3 . The 
velocity-averaged annihilation cross section is taken as 
(cr aim w) = 3 x 10~ 26 cm 3 s _1 , with an annihilation chan- 



nel into b-b. According to 3l|, the propagated primary 
p flux is only very mildly dependent on the annihilation 
channel and the DM distribution function. Therefore, 
our assumptions can be considered valid for a generic 
WIMP dark matter candidate except for a rough rescal- 
ing factor. Propagation is treated in the same way as 
for the secondary component 29j, |3l|. As a reference 
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FIG. 2: Upper limits on the enhancement factor to the pri- 
mary p flux as a function of the WIMP mass, derived from 
a comparison with PAMELA high energy data. Each curve 
is labelled according to the corresponding PAMELA energy 
bin. 



case, we employ the best fit transport parameters listed 
above and recall that the uncertainty on the primary p 
flux due to propagation spans roughly one order of mag- 
nitude above and one below the best fit scenario. We add 
the calculated primary p flux for different WIMP masses 
to the secondary component and compare the total flux 
to PAMELA high energy data, namely T p > 10 GeV. To 
be conservative, the background calculated from Bring- 
mann & Salati's p production cross sections is considered 
(dashed curves in Fig. QJ. We derive the factor by which 
the DM flux could be enhanced without exceeding ex- 
perimental data (2a error bars) in any energy bin. The 
maximum allowed enhancement factor is plotted in Fig. [2] 
as a function of the WIMP mass: it cannot exceed 6-20- 
40 for towimp=100-500-1000 GeV, respectively. These 
limits can be reinforced as well as relaxed by quite simple 
modifications of the key ingredients in the flux calcula- 
tion, just as described above. The boost factor may be 
ascribed, in principle, to dumpiness in the DM distribu- 
tion [34I ] — this contribution being energy-dependent — as 
well as to an increase of the annihilation cross section 



as proposed by 35[ and more recently by [36j using the 
Sommerfeld effect. 

Our conclusions have important consequences on the 
explanations of the positron data based on the annihi- 
lation of DM species within the Milky Way halo. The 
positron fraction suffers from large uncertainties related 
for instance to the poorly determined electron spectral 
index above 10 GeV Although soft electrons are as- 



sociated to large values of the positron fraction and to a 
marginal agreement of the pure secondary positron flux 
with the measurements, wc cannot dismiss the possibil- 
ity of a hard cosmic ray electron distribution. A spectral 
index of 3.44 4l| leads actually in the left panel of Fig. [3] 
to the long-dashed curve featuring a low background 
case. With a typical annihilation cross section (cr a nnw} 
of 3 x 10~ 26 cm 3 s _1 , WIMPs do not produce enough 
positrons to reproduce the increasing trend observed in 
e + 1 (e + + e~) data [401 ] ,so that a significant enhancement 
of the annihilation rate is necessary as shown in [iil ]. 
However the boost factor associated to DM clumps can- 
not exceed at most a factor of ~ 10 in the standard A- 
CDM scenario of structure formation [34|. Astrophysics 



does not provide then a natural explanation for the large 
boost factors required to fit the positron excess. That is 
why the Sommerfeld effect [HB] has been advocated as a 
plausible mechanism to significantly increase the WIMP 
annihilation cross section in the non-relativistic regime 
prevailing today in galactic haloes. Heavy DM species is 
a prerequisite. We then consider a generic 1 TeV parti- 
cle annihilating into W + W~ pairs and boost (a ann v) by 
a factor of 400 in order to get the solid line in the left 
panel of Fig.[3l Although an annihilation cross section of 
1.2 x 10 -23 cm 3 s _1 is possible should non-perturbative 
effects be involved, the consequences on antiprotons are 
drastic. The red solid curve in the right panel of Fig. [3] 
features an unacceptable distortion of the p spectrum. 
The DM positron signal cannot be enhanced without 
playing havoc with the p measurements. 

Nonetheless, notice that ways out are possible whose 
careful investigation is beyond the scope of this Let- 
ter. The value of 400 assumed for the positron signal 
of Fig. [3] could arise from the combined effects of DM 
dumpiness and (<T ann w) enhancement. If a generous fac- 
tor of 10 is assumed for the former — a marginally accept- 
able value — the latter does not exceed 40. Unlike 
positrons which are produced locally, the antiprotons de- 
tected at the Earth originate from a large region of the 
Milky Way halo over which substructures may not be as 
important as in our vicinity. The p flux may not be much 
enhanced by the presence of DM clumps so that a value of 
40 would apply in that case to the antiproton boost. The 
corresponding blue long-dashed line in the right panel of 
Fig. [3] features a fairly acceptable p spectrum. 

Viable scenarios such as a large black hole population 
pervading the Galaxy [43j, |44[ also lead to large boost 
factors although it seems difficult a priori to prevent the 
antiproton flux from being enhanced too much. Notice 
finally that the cosmic ray propagation model could be 
different from the one selected here. Once again, positron 



4 




Modulation with F = 500 MV 
Scan with B/C compatible data 



1 10 10 2 10 o.l 1 10 100 1000 

Energy [GeV] tj» [Gev] 



FIG. 3: The fiducial case of a 1 TeV LSP annihilating into a W + W pair is featured. In the left panel, the positron signal 
which this DM species yields has been increased by a factor of 400, hence the solid curve and a marginal agreement with the 
PAMELA data. Positron fraction data are from HEAT AMS-01 [HIH and PAMELA 0. If the so-called Sommerfeld 
effect is invoked to explain such a large enhancement of the annihilation cross section, the same boost applies to antiprotons 
and leads to an unacceptable distortion of their spectrum as indicated by the red solid line of the right panel. 



and antiproton fluxes have different behaviors toward a 
change in the propagation parameters. For example, the 
primary p flux could be easily decreased by an order of 
magnitude without violating B/C data, allowing a Som- 
merfeld boost of the cross section of 400. 

A new calculation for the secondary cosmic antiproton 
flux and the relevant uncertainties have been presented. 
The ratio p/p has been derived after fitting recent pro- 
ton data. Our predictions can explain the experimental 
data, and in particular the recent PAMELA data, which 
span more than two decades in energy. No exotic contri- 
bution — as from annihilating dark matter in the galactic 
halo — has to be invoked to reproduce experimental re- 
sults. Analyzing the high energy part of the PAMELA 
p/p we derive strong upper limits on possible enhance- 
ments of the exotic p flux as a function of the WIMP 
mass. Relying on standard assumptions the exotic an- 
tiproton flux induced by a to wimp = 100 GeV (1 TeV) 
DM halo cannot be increased by more than a factor 6 (40) 
without overrunning PAMELA data. Would the Som- 
merfeld effect (PF+PF - channel) be invoked to explain 
PAMELA leptonic data, the corresponding enhancement 
of the p production would lead to an unacceptable dis- 
tortion of the p/p spectrum. 
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